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A method is presented that, when used in conjunction with single molecule experimental tech- 
niques, allows for the extraction of rates and mechanical properties of a biomolecule undergoing 
transitions between mechanically distinct states. This analysis enables the exploration of systems 
where the lifetimes of survival are of order of the intrinsic time constant of the experimental ap- 
paratus; permitting the study of kinetic events whose transition rates are an order of magnitude 
(or two) larger than those that can be studied with traditional averaging methods. Using cur- 
rent experimental techniques, this allows for the study of biomolecules whose lifetime of survival 
in a given state are as low as milliseconds down to microseconds. The relevant observable is the 
auto-correlation function of the experimental probe that is attached to the biomolecule of interest. 
Obtaining these auto-correlations functions is challenging and our analysis is based upon a series 
solution. Although the general solution is expressed in terms of a series, two important limits are 
studied, whose solutions are obtained exactly. These limits correspond to physically opposing lim- 
its: where transitions between the mechanically distinct states of the biomolecule occur on either 
much faster or much slower time scales than those governing the motion of the experimental probe. 
Not only do these limits correspond to physically opposing bounds, but it is also proven that they 
correspond to mathematical bounds to the general solution. Motivated by the derivation of these 
opposing bounds, two series solutions for the general case are then presented. Armed with these 
two general series solutions that approach the exact solution from opposing points, we then present 
an error analysis for truncating each series at arbitrary order and obtain a range of parameters 
for which this method could be applicable to the study of the two state biomolecular problem. 
Finally, both series (up to third order) are expressed for the two state problem when the system 
obeys Markov statistics. These solutions should be amiable to the analysis of experimental data, 
expanding temporal analysis of data from single molecule experiments. 

PACS numbers: 02.50.Fz 82.37.Np 87.15.La 82.37.Gk 



I. INTRODUCTION 

The advent of experimental techniques such as atomic 
force microscopy (AFM) and optical- and magnetic- 
trapping has made it possible to manipulate individual 
biological moleculesjl B 

El El El mill Ell ■ These methods provide comple- 
mentary insights to those obtained using the techniques 
of solution biochemistry. In particular, they allow for di- 
rect interpretation of specific processes occurring on in- 
dividual molecules which include the fluctuational char- 
acter of real macromolecular trajectories. 

The single molecule experiments that have measured 
the kinetics of structural transitions of macromolecules 
have accessed a wide variety of processes. Some of the 
most compelling examples include the determination of 
rates, including those associated with the folding and un- 
folding of various proteins HUE! El El and RNA [H, 
structural transitions of polysaccharides p^ . protein- 
induced loop formation in DNA 0, ^| and the binding 
and unbinding of myosin to actin filaments 0, H IT^ [l8l | . 
Despite the diversity of these experiments, many have 
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been limited to the study of kinetics where the life- 
times of the macromolecular states of interest are much 
larger than the intrinsic time constant of the experimen- 
tal probe. 

As an illustrative example of this aforementioned lim- 
itation and to exemplify the types of problems that in- 
spired this work we now consider the abstract case where 
a macromolecule is attached to a force probe (AFM, for 
example) and alternates between two mechanically dis- 
tinct states. The inset in Fig. ^ displays the problem of 
interest, where the macromolecule (DNA, for example) 
of concern interacts with a protein that binds to specific 
sites (diamonds in inset) on the macromolecule and loops 
out the intervening fragment. The macromolecule then 
fiuctuates between a looped and unlooped state. (One 
can also consider folded/unfolded states or other types of 
macromolecular interactions which cause conformational 
changes.) From the point of view of the experiment, 
these states can be distinguished only if they give rise 
to different probability distributions for the observable 
associated with the probe. An example of this scenario 
is displayed in the figure where the probability of finding 
the probe displaced by a certain amount (Pa or Pb : thick 
solid lines in the figure) depends upon the state of the 
macromolecule (state A or B). (In the figure it has been 
assumed that the distribution functions are Gaussian.) 
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FIG. 1: Inset: System of inspiration. An AFM is attached 
to a macromolecule which can exist in two states, an un- 
looped state (A) or a looped state (B). The looped state is 
caused by an interaction between the macromolecule and a 
protein which binds to specific sites (diamonds) on the macro- 
molecule. Figure: Different probability distribution functions 
(P(a;)) of the displacement of the AFM tip exists for each 
state (Pa (a;), PB(a;): thick solid lines). For cases of interest 
the instantaneous distributions (P(2;)) overlap. If the position 
of the probe is averaged ({x)t) over some time window T, 
then the probability distribution function (P((x)t): dashed 
lines) of this averaged value will separate, provided T is large 
enough and the likelihood of transitions during T is small. 
Pa,b (thin solid line): instantaneous joint probability distri- 
bution function of the system. 



Deleterious to the problem is the fact that in order to 
have transitions between the macromolecular states the 
probe must access positions where the molecule can ex- 
ist in cither state. Therefore, the distribution functions 
must overlap and in many cases it is imperative that the 
distribution functions overlap sufficiently in order that 
numerous transitions occur during an experimental run. 
Because of the requirement for overlap a single measure- 
ment, in general, cannot determine the state of the sys- 
tem. To overcome this dilemma measurements x (x: dis- 
placement of the probe for example) are averaged (x)t 
over a time window T which reduces thermal noise. If T 
is sufficiently large the probability distribution functions 
of the averaged value (P((x)t): dash-lines in figure) will 
separate and the averaged measurements ((a;)T) can dis- 
cern the states of the system. (Alternatively, moments 
other than the first {{x)t) can be used to discern the 
states of the system. The arguments that follow here 
apply to such measurements (see Sect. 111)1.) 

Given this framework we now illustrate the aforemen- 
tioned limitation that exists within single molecule ex- 
perimental studies which, when using such averaging pro- 
cedures, requires the lifetimes of the states to be much 
larger than the intrinsic time constant of the experimen- 



tal probe. This limitation results from the fact that 
in order to obtain well separated distribution functions 
(P((x)t)) which properly correspond to the states of the 
system requires a relationship between the time window 
T of averaging and the two aforementioned time scales. 
First, in order that the resulting distribution functions 
P((x)t) of the averaged value {x)t correspond to a par- 
ticular state requires that the likelihood of a transition 
from one state to the other is very low over the time 
period T. This is guaranteed only if the lifetimes of the 
states {Ti} (average survival time) is much greater than 
the measuring period T, {%} ^ T. Secondly, in order 
to obtain well separated distributions requires that suf- 
ficient thermal noise is reduced over the time window T. 
This is satisfied if the intrinsic time constant of the exper- 
imental probe Tp is much smaller than the measuring pe- 
riod T, Tp <C T. (For an over-damped system rp = ^/kp, 
where 7 is the damping constant of the probe and fcp is 
its spring constant.) The argument for this relation goes 
as follows: Given an instantaneous measurement of the 
probe Xi which lies in the overlap region, for the aver- 
aged value {x)'T (for which this instantaneous measure- 
ments Xi contributes) to fall in one of two well separated 
distribution functions requires the set of instantaneous 
measurements {a;} (which determines this averaged value 
(x)t) to represent a distinguishing characteristic of that 
particular distribution function. This requires the set of 
instantaneous measurements {x} to contain a range of 
values, implying it must contain a number of uncorrc- 
lated measurements. Roughly speaking, it takes time of 
order of the intrinsic time constant of the experimental 
probe Tp to make two uncorrelated measurement; there- 
fore, requiring rp <C T. Figure ^ displays the distribution 
functions of the averaged value, P((x)t), corresponding 
to a time window with ten uncorrelated measurements, T 
« 10 Tp. Under such conditions two "separated" distri- 
butions functions are obtained. The quantitative reduc- 
tion of the width with respect to the number of uncor- 
related measurements is presented in Appendix ^ (In 
principle, the time constant governing the motion of the 
probe should include contributions from the molecule, 
however, for the systems of interest the intrinsic time 
constant of the experimental probe Tp gives the appro- 
priate measure. This point is further elaborated on in 
Appendix ^ 

These requirements have placed limitations on nu- 
merous experimental studies that use single-molecule 
probe techniques: the lifetimes of survival must be 
much larger than the intrinsic time constant of the ex- 
perimental probe |7i| ^ Tp. Optical trap measure- 
ments 0, nil, 111 111 111 El typically have a time con- 
stant of the order of 0.5 ms. Such techniques have been 
limited to the study of kinetic events where states survive 
on times scales of the order of lO's of milliseconds and, 
in general, much longer. H H H E E El 

A number of biological molecules can fluctuate between 
different states on time scales of the order of microsec- 
onds down to nanoseconds. Such measurements have 
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been made usin g th e techniques of solution biochem- 
istry [23, 1^ 123. l23l| . For a single- molecule probe tech- 
nique to have the capability of measuring the properties 
of such fast switching systems requires the development 
of either new experimental methods with higher tempo- 
ral resolution or alternative analytic methods capable of 
extracting properties of the system from current experi- 
mental measurements. New experimental methods such 
as the placement of nanometer sized resonators on stan- 
dard AFM tips [23. |25| can provide a promising future 
in the development of probes with both the high tempo- 
ral resolution and the large compliance needed to mea- 
sure such properties. However, such systems have not yet 
been fully developed for the study of single biomolecular 
systems. 

In this work we propose an alternate analytic approach 
to analyze data from single- molecular probe techniques. 
This analysis allows for the exploration of macromolecu- 
lar systems when the lifetimes of survival of the states are 
on the order of the intrinsic time constant of the experi- 
mental probe, {%} ~ Tp. The property computed is the 
position-position auto-correlation function of the exper- 
imental probe. The auto-correlation function is chosen 
because it depends upon the joint probability distribu- 
tion Pa.b (thin line in Fig. ^| function and not on the 
distribution function of the individual states (Pa or Pb). 
This allows us to forgo the requirement that the system 
remain in a particular state for a "long" period of time 
and in turn increases the time resolution of the analy- 
sis. (The trade-off, however, is in the ability to relate the 
statistical properties of the entire system (one which de- 
pends upon Pa,b) to properties of the individual states 
(lifetimes for example).) Moreover, the auto-correlation 
function varies on time scales of order of the intrinsic 
time constant of the experimental probe; therefore, its 
form should be affected by events (macromolecular tran- 
sitions) that occur on such corresponding time scales. By 
appropriate analysis of the data, it should be possible 
to determine lifetimes of states with characteristic scales 
ranging from milliseconds to microseconds with current 
experimental techniques. 

This manuscript focuses on the general mathematical 
formalism for calculating the auto-correlation function 
of a system that is switching between mechanically dis- 
tinct states. To further illustrate the problems of in- 
terest we proceed in Sect. ^ with numerical studies of 
a mock experiment (computer simulation) representing 
the systems of interest. These simulations reveal how 
traditional methods (such as that described in this in- 
troduction) fail to determine the rates of transitions of 
a biomolecule when fluctuations occur on times scales of 
order of the intrinsic time constant of the experimental 
probe. In addition, presented arc numerical solutions for 
the auto-correlation function of the mock system. These 
numerical solutions illustrate how its form depends on 
the value of the rates of transition and could be used as 
a conduit to extract such rates from experimental mea- 
surements, those beyond the reach of traditional analysis. 



The following sections proceed with the mathemat- 
ical formalism associated with calculating the auto- 
correlation function. Prior to determining the general 
solution we focus on two special cases, termed the fast- 
switching and slow-switching limits. Sect. IIIII These 
limits correspond to cases in which the macromolecule 
switches between states on time scales that are either 
much faster or much slower than the time scales associ- 
ated with the motion of the experimental probe. In ad- 
dition, "closed" form solutions for auto-correlation 
function in these limiting scenarios are presented and 
it is proven that these solutions correspond to either a 
lower bound (fast-switching) or an upper bound (slow- 
switching) to the general solution for the auto-correlation 
function. 

Although "closed" form solutions are found for these 
two special cases, they do not exists in general, at least 
to our knowledge. Therefore, solutions are provided for 
the auto-correlation function in the form of a series solu- 
tion. Sect. II VI Motivated by the fact that the two bounds 
found for the general solution also correspond to physi- 
cally opposing limits, we generate two series solutions for 
the auto-correlation function which initiated at these op- 
posing bounds. Providing these two series solutions has 
the appeal that they approach the exact solution from 
opposite points. 

The convergence of these two series to the exact solu- 
tion is examined in Sect.0 Additionally, this section ad- 
dresses issues regarding the applicability of this method 
to the analysis of experimental data. This section focuses 
mainly on the two state problem and provides solutions 
for both series, truncated at third order, for the two state 
Markov process. These solutions should be amiable to ex- 
perimental analysis, where the rates of transition can be 
determined in a regime inaccessible by traditional analy- 
sis. 



II. ILLUSTRATIVE CASE: A NUMERICAL 
STUDY 

Our work is motivated primarily by biological 
molecules which can make transitions between mechan- 
ically distinct states and the associated single-molecule 
probe techniques used for detecting the properties of 
these systems. In general, the transitions we have in mind 
correspond to structural transitions of the biomolecule 
which are induced either by thermal fluctuations or im- 
posed by interactions with other macromolecules. (See 
inset in Fig. ^ for example.) Figure [21 presents a 
schematic of the energy landscape describing the con- 
formational states of a biomolecule which can exist in 
any of a number of different folded states. (The energy 
landscape may be thought of as the bare energy land- 
scape of the molecule or the resulting energy landscape 
which includes to the presence of the experimental probe 
and is generally biased along a particular reaction co- 
ordinate due to the stretching of the molecule.) Each 
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Reaction Coordinate 

FIG. 2: Energy landscape of a biomolecule undergoing struc- 
tural transformations between folded states. The reaction 
coordinate corresponds to a physical degree of freedom along 
which the biomolecule can make a transition from one state to 
the others. Each folded state corresponds to a well in the en- 
ergy landscape and their associated mechanical response are 
distinct. Because of thermal fluctuations the biomolecule can 
make transitions between the different states. 



folded state has a different mechanical response and cor- 
respond to distinct energy wells in the energy landscape 
of the biomolecule. Because the biomolecule is undergo- 
ing thermal fluctuations it can make transition between 
the various local wells. 

In single molecule experiments like those being consid- 
ered here, an experimental probe is attached to the end 
of a biomolecule. Each detectable state differs in their 
mechanical response. The transitions between the states 
are governed by some set of rate constants Kij , where the 
subscripts label the state being exited and that being en- 
tered. Figure 13 (left) displays a typical example. The 
experimental probe (optical bead in an optical trap, for 
example) is held essentially at a fixed position (other than 
motion induced by thermal noise) and the biomolecule is 
essentially held at fixed extension. The mechanical char- 
acteristics of the biomolecular states will reveal them- 
selves as mechanical springs whose corresponding spring 
constants are determined by the state of the system at 
that given extension (2^. The plot on the right displays 
a cartoon representing the mapping of the original exper- 
imental scenario onto a corresponding mechanical model 
in terms of coupled masses and springs. The experimen- 
tal apparatus and the molecule in its different states have 
mechanical stiffnesses characterized by spring constants 
kp and fcmi, respectively, where fc,„j; differs for each molec- 
ular state i. 

As expressed in the introduction, the states of the sys- 
tem can be determined, in general, only if the lifetimes 
of the states are sufficiently large as compared to that 
of the intrinsic time constant of the experimental probe, 
{Ti} ^ Tp. To give a more concrete representation of 
this limitation we now give a specific study by way of nu- 
merical methods, where we perform mock experiments 
to determine the feasibility of resolving the properties of 
the different states of a macromolecule through meth- 
ods similar to that presented in the introduction. Next, 






FIG. 3; An experimental scenario. Left: An optical bead 
in an optical trap is connected to a biomolecule which can 
access multiple, mechanically distinct states. The other end 
of the biomolecule is anchored in place. Transitions between 
the states are determined by a set of rate constants {^ij)- 
Right: Schematic of interactions in terms of springs. The 
probe (optical bead) has an intrinsic spring constant (fcp) and 
each biomolecular state has a speciflc spring constant (fcmi)- 
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FIG. 4: Division between the over-damped and under- 
damped regimes as a function of frequency lo and viscous 
damping 7/2M, where 7 is the damping constant and M is 
the associated mass of the probe. Region where optical trap 
(OT) experiments lie is encased in the dashed-circle. The 
dashed line is the separation between the over-damped and 
under-damped regions. 



we present numerical solutions for the auto-correlation 
function to demonstrate how an analysis of the auto- 
correlation function can be used to alleviate this limi- 
tation. Having demonstrated the utility of analyzing the 
auto-correlation we proceed this section with our theoret- 
ical treatment for calculating the auto-correlation func- 
tion. 

Prior to proceeding with this illustration expressed are 
a few key assumptions imposed on the system, which 
in turn affects the form of the equation of motion em- 
ployed for the experimental probe. These assumptions 
are physically motivated and will be used throughout the 
remainder of the paper. First, it is assumed that the 
different states of the biomolecule only affect the spring 
constant of the biomolecule and do not alter the damp- 
ing constant of the experimental probe. This assump- 
tion is justified on the grounds that the experimental 
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probe is, in general, much larger than the molecule of 
interest and consequently, any fluid flow generated by 
the biomolecule (or more importantly the differences in 
fluid flow from different states) will only weakly influ- 
ence the experimental probe. As a result, such differences 
are ignored. Additionally, it is assumed that the experi- 
mental probe is over-damped. For most single molecule 
experiments of interest this will be the case. Figure ^ 
charts the dividing line between the over-damped and 
under-damped regimes for optical trap measurements, as 
a function of frequency of the probe (w) and its viscous 
damping (7/2M), where 7 is the damping constant asso- 
ciated with the motion of the probe and M is its mass. 
Clearly these measurements are performed in the over- 
damped limit. Furthermore, it is assumed that once a 
state (in one energy well) begins a transition to another 
state (to another energy well) it occurs instantaneously 
or more specifically, faster than all other time scales in 
the problem; allowing the mechanical properties of the 
macromolecule at any given instant to be described by a 
state corresponding to one of the energy wells. In addi- 
tion, it is assumed that the interaction between the heat 
bath and experimental probe can be expressed as a white, 
Gaussian noise process. Finally, in this manuscript it is 
assumed that ergodicity holds, time averages are equal 
to ensemble averages. 

A. Limitations in averaging procedures 

We now present our numerical illustration. Given the 
above assumptions, the equation of motion for the probe 
is written using a Langevin approach, in which the force 
acting on the probe exhibits switching behavior charac- 
terized by different spring constants. It is mathematically 
convenient to think of this problem as a single Langevin 
equation with a constant damping constant, but with a 
fluctuating spring constant and equilibrium position. Ex- 
pressed is the equation of motion for the tethered probe: 

7i(<) = Frit) - k{t){x{t) ~ x°{t)). (1) 

Here x{t) is the position of the probe and the associated 
velocity is given by x{t). The fluctuating spring con- 
stant is represented by k{t), whose value takes ki when 
the system is in state i. The value for the fluctuating 
spring constant, ki, includes both the spring constant of 
the biomolecule and the intrinsic spring constant of the 
experimental probe 

ki = fcp + A^mz- (2) 

The equilibrium position of the system x°{t) is also a 
fluctuating variable, whose value is denoted as x° when 
the system is in state i. The damping constant, 7, is 
the same for each state [s^. Finally, Fr{t) is the random 
thermal force due to the presence of the environment and 
is assumed to be Gaussian white noise characterized by 



the correlations 

{Fr{i)) = 0, (3) 

{Fr{Q)Fr{i)) = 2k^T-i5{t), (4) 

where ke is the Boltzmann constant and T is the tem- 
perature. 

To develop intuition for this problem, we choose 
the spring constant for the experimental probe fcp to 
be 20 pN//im and its damping constant 7 to be 9.5 
fg/ns. Such values arc typical of optical beads in a 
laser trap jl^ . The corresponding intrinsic time con- 
stant (rp = l/k-p) of the experimental probe equals 0.48 
milliseconds. We consider a biomolecule which can ex- 
ist in two states with corresponding spring constants 
25 pN/^m (state 1) and 5 pN//xm (state 2). Spring 
constants which are easily obtainable in muscle pro- 
teins 0, 0, O, ■ For simplicity, the equilibrium 
position of each state is chosen to be equal, corresponding 
to observation of transverse motion of the experimental 
probe. (Transverse to the extension of the biomolecule.) 
The rate for the biomolecule to make a transition from 
state i to state j is given by the rate constant . 

Before embarking on our numerical test we investi- 
gate properties of the distribution function for the po- 
sition of the experimental probe when the molecule is 
in a given state. These distribution functions (pi{x)) arc 
determined through Boltzmann statistics, each with form 

Pii.x) . 

Here /3 — l/keT is proportional to the inverse temper- 
ature and Zi is the partition function (normalization) 
for state i. (Recall ki is the spring constant for state i, 
Eq-El) Figure |31 displays the two distribution functions. 
Unlike the example in the introduction, here the average 
position of the probe ^x)i (first moment) for each state 
cannot be used to decipher the state of the system as 
they are equal, {x)\ ~ {x)2- In such cases the squared 
displacement of the probe {x^)i (second moment) differ 
for each state (due to differences in spring constants) and 
can be used to decipher the state of the system. The 
squared displacement corresponds to a measurement of 
the width of a distribution function squared. Such pro- 
cedures have been employed in single molecule experi- 
ments nil El El 113. 

Although the measure used here differs from that pro- 
posed in the introduction, the arguments from the in- 
troduction arc still applicable in revealing constraints 
on various time scales. To obtain an accurate measure- 
ment for the width of a distribution function requires the 
time window of averaging T to be sufficient, such that a 
number of uncorrelated measurements are made within 
this time window; requiring the time window to be much 
greater than the intrinsic time constant of the experimen- 
tal probe, T3> Tp. To obtain accurate accounts for the 
individual widths of the distribution functions requires 
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position X [nm] 

FIG. 5: The normalized distribution functions of the posi- 
tion of the experimental probe when the molecule is in ei- 
ther state 1 or state 2. The spring constant of each state 
{k = kp + kmi) is equal to the spring constant of the exper- 
imental probe (fcp) plus the spring constant of the molecule 
for that particular state (fcmi). For state 1 the molecular 
spring constant is 25 pN//xm and for state 2 it corresponds to 
5 pN//im. The spring constant of the experimental probe is 
taken to be 20 pN//im. 



TABLE I: Rates of transition between the two different 
states. Kij is the rate at which transitions are made from 
state i to state j. The rates correspond to a range of illus- 
trative cases, which differs by powers of ten. The rates are in 
units of Hertz and in the parenthesis are the lifetimes of the 
states (inverse of the rates) in units of milliseconds. 



rates 


case A 


case B 


case C 


case D 


Ki2 [Hz] (ms) 
K21 [Hz] (ms) 


0.26 (3800) 
0.11 (9100) 


2.6 (380) 
1.1 (910) 


26 (38) 
11 (91) 


2600 (0.38) 
1100 (0.91) 



the system to remain in a given state over this time win- 
dow; requiring the time window to be much less than the 
Ufetimes of the states, {^}- These two conditions 
lead to the requirement that for this accounting method 
to be useful for exploring the states of the system the 
intrinsic time constant of the experimental probe must 
be much less than the lifetimes of the states, rp ^ {%}■ 
Having illustrated that the time constraints presented 
in the introduction still pertain to our numerical study, 
we return to our numerical test and illustrate the limi- 
tations of a traditional analysis for identifying the states 
of the macromolecule and determining their associated 
properties. To explore the applicability of such tradi- 
tional methods we considered a range of examples. These 
examples range with respect to lifetime of the states, from 
those which are much greater than the intrinsic time con- 
stant of the experimental probe {{%} ^ rp = 0.48 ms) 
to those which are of the same order as the intrinsic time 
constant of the experimental probe {{%} w rp). Four 



cases are considered within this range, whose rates are 
displayed in Table The lifetimes of the states are dis- 
played in parentheses in the table and cover four orders 
of magnitude. 

To obtain a numerical solution for the position of the 
probe, x{t), Eq.^is i nteg rated using a standard Langevin 
dynamics algorithm pa | [sol j . During the calculation the 
state of the system and hence the spring constant (fc(i)) 
is determined as follows; If at time t the system is in 
state i then the probability of the system at time t + dt 
to be in state j is Kijdt. Transitions are then determined 
by a uniform random number generator, which generates 
a number < r < 1. If r < Hijdt then a transition 
occurs. Such an algorithm corresponds to a two state 
Markov process and gives the proper statistics so long as 
dt is chosen such that Kijdt <^ 1, for i ^ j. 

Given the above algorithm, the equation of motion is 
integrated and the position of the probe is determined. 
Initially case A is only considered, where the lifetime of 
survival is much greater than the intrinsic time constant 
of the experimental probe. Here, traditional methods 
(like those discussed here) are applicable. 

Because instantaneous measurements cannot decipher 
the state of the system, identification of the states have 
relied on averaging procedures. To illustrate the need 
for some type of averaging procedure the instantaneous 
squared displacement of the probe x^t)"^ is displayed in 
Fig. 13 (gray dots). Instantaneous measurements do not 
provide enough information to reveal the state of the sys- 
tem and are not accurate predictors for the mean squared 
displacement of the probe. The mean squared displace- 
ment for state i corresponds to {x'^)i = kBT/(fcp -I- k„^i) 
and is displayed as a black solid line in the figure. 

As argued, when averaged over a number of uncorre- 
latcd measurements M (M w T/rp) can accurate pre- 
dictions for the mean squared displacement (second mo- 
ment) be obtained, revealing properties of the biomolec- 
ular state. In addition, if this time interval small com- 
pared to the lifetimes of the states (rp ^ T ^ {^}) then 
the individual states of the system can be identified. To 
exemplify this point Fig. [T] (a) displays the squared dis- 
placement of the probe {x{t)'^)T, for case A, where the 
squared displacements are averaged over a 50 millisecond 
time interval (number of uncorrelated measurements M 
~ 100). Here the gray dots correspond to the time aver- 
aged value (x(t)^)T while the black solid line corresponds 
to the mean squared displacement {x'^)i for that particu- 
lar state. Unlike the results in Fig. El now the state of the 
molecule can be identified at time t and the properties of 
the states can be determined. The mechanical properties 
can be determined by calculating the squared displace- 
ment in each state and setting it to kBT/(fcp -|- fcmi) and 
the lifetimes of the states can be determined by accumu- 
lating the survival time of each state. 

We now explore cases where the rates are increased by 
orders of magnitude in order to determine the limitations 
of the above averaging procedure. Initially, the rates are 
increased by an order of magnitude, corresponding to life- 



FIG. 6: The figure displays the instantaneous squared dis- 
placement {x{t)^) of the experimental probe, gray dots. The 
rates of transition between the states are K12 = 0.26 Hz and 
K21 ~ 0.11 Hz. The solid bold line displays the mean squared 
displacement {{x'^) = kBT/(fcp + kmi)) of the probe for the 
given state. 



FIG. 7; The squared displacements averaged, {x{t)^)T (gray 
dots), over a small time interval T. The solid bold lines corre- 
spond to the mean squared displacement {{x'^) = kBT/(fcp -|- 
fcmi)) of the probe for the given state, (a) Case A: /ti2 = 
0.26 Hz, K21 = 0.11 Hz and T = 100 milliseconds, (b) Case 
B: K12 = 2.6 Hz, K21 ~ 1.1 Hz and T = 50 milliseconds, (c) 
Case C: K12 = 26 Hz, K21 ~ 11 Hz and T = 5 milliseconds, 
(d) Case D: K12 = 2600 Hz, K.21 = 1100 Hz and T = 0.05 
milliseconds. 



times of 380 ms and 910 ms (case B). Figure [7|(b) dis- 
plays the averaged squared displacement (x(t)^)T where 
the time window for averaging T is taken to be 30 ms 
(M « 60). Here too, we are in the regime where the life- 
times of the states are much greater than the intrinsic 
time constant of the experimental probe and; therefore, 
a time window can be taken such that the states can be 
accurately identified and the properties of the states can 
be determined. Results for increasing the rate by an- 
other order of magnitude, corresponding to lifetimes of 
38 ms and 91 ms (case C), are displayed in Fig. [3 (c). 
The squared displacement is averaged over a time win- 
dow T = 5 ms (M w 10) which is an order of magnitude 
larger than the intrinsic time constant of the experimen- 
tal probe. For this case it becomes extremely difficult to 
identify the states of the system. Finally, case D is dis- 
played in Fig. (d), where the rates are now increased 
by an additional two orders of magnitude. Here the rates 
are of order of the inverse intrinsic time constant of the 
experimental probe. The squared displacements are av- 
eraged over a time window T = 0.05 ms (M « 1). Here 
too, the states of the system cannot be resolved when 
using this conventional technique. 

The failure to identify the states and their properties 
in case C (Fig. [Jc)) and case D (Fig. Cfd)) lies in the 
fact that the time window of averaging T is not large 
enough to contain a significant number of uncorrelated 
measurements. It may be envisioned that this failure 
could be alleviated if the time window of averaging T 
was increased. However, to accurately identify the states 
of the system (and hence determine their mechanical and 
kinetic properties) the likelihood of a transition occurring 



within this measuring period T must be small. This re- 
quires a time window an order of magnitude less than 
(and generally much less than) the lifetime of the states, 
T <C {%}■ Increasing the time window will result in a 
number of measurements corresponding to mixed states, 
not providing better statistics regarding the properties of 
the individual states. 

Another possible route in attempting to alleviate the 
failure, in identifying the properties of the states for case 
C and D, might rely on increasing the frame rate of obser- 
vation. By frame rate of observation we mean the rate at 
which the instantaneous measurements are stored. Upon 
increasing results in a larger set of data points to calcu- 
late the second moment (a;(t)^)T over fixed time window 
T. However, increasing the frame rate of observation 
docs not reduce the noise, because it does not increase the 
number of uncorrelated measurements. Many uncorre- 
lated measurements are needed to obtain good statistics 
and it takes time of order Tp to make two uncorrelated 
measurements. When averaged over the same time win- 
dow T the features of the plots in Fig.|7|are not altered 
when the frame rate of observation is increased. 

We have illustrated that to resolve the dynamics and 
statistics of a biomolecule, which is undergoing structural 
transitions, requires an experimental apparatus whose in- 
trinsic time constant is much smaller than the lifetimes 
of the states, rp <^ {%}■ This requirement is based 
upon the precondition that the properties of the indi- 
vidual states are to be identified separately; requiring an 
experimental probe to have sufficient temporal resolution 
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to adequately explore each state prior to the occurrence 
of a transition. One of the key ambitions of the remain- 
der of the paper is to develop a formalism that allows 
us to forgo this precondition and in turn alleviate the 
demand that Tp ^ {%}■ 

To alleviate this demand we propose an alternative 
analytic approach for analyzing the data from single- 
molecule probe techniques. The property that we com- 
pute is the position-position auto-correlation function. 
To illustrate the potential benefit of analyzing the auto- 
correlation function, we now present numerical solutions 
for sample cases presented in this section. This analysis 
will lend insight into how the form of the auto-correlation 
function is influenced by the rates of transition between 
biomolecular states and how such properties are revealed. 
Following this analysis we present our general mathemat- 
ical theory. 



B. Numerical solution for the auto-correlation 
function 

We present numerical solutions for the auto-correlation 
function for two of the cases just analyzed, case B and D. 
The auto-correlation function is defined as 



{Sx{t)Sx{0)) = {xit)x{0)) - {xf 



(5) 



where 5x{t) is the displacement of the probe from the 
equilibrium value {{x)) at time t, 

dx{t) = x{t) - (x). 

The statistical averages in Eq.jS] can be interpreted in 
terms of time averages. The first term on the right-hand 
side is defined as 

{x{t)x{0)) = hm i r dt' x{t + t')x{t') (6) 

T — OO i Jg 

and the equilibrium position of the probe is defined as 



(^) = 1™ T / dt' x{t'). 



(7) 



The auto-correlation function measures the memory of 
particle displacements. Given an initial displacement at 
time t = 0, the auto-correlation function gives the time 
scale at which the particle returns to its equilibrium value 
and looses accounts of that given initial displacement. If 
the likelihood of a transition between biomolecular states 
is high during this relaxation period, then upon return 
to its equilibrium position the particle itself should have 
memory regarding the various states visited during its 
journey; therefore, such transitions should reveal them- 
selves through the solution to the auto-correlation func- 
tion. The auto-correlation function decays on time scales 
of the order of the intrinsic time constant of the experi- 
mental probe, potentially capable of revealing transitions 




time [ms] 



FIG. 8: Auto-correlation function for a system with a fluc- 
tuating spring constant. Case B and case D are represented 
by blue circles and red circles, respectively. The solid lines 
refer to physically motivated guesses for the solution of the 
auto-correlation function. These guesses are expressed in 
Eqs.lHland^l blue line and red line, respectively. 



occurring on these time scales. 

To explore this possibility, we now present numerical 
solutions for the auto-correlation function for case B and 
case D. Given the numerical solution for the position of 
the probe {x{t)) the auto-correlation function can be de- 
termined through application of Eq.|31 The results are 
displayed in Fig. |S1 where the solution for case B is dis- 
played as blue circles and that for case D is displayed 
as red circles. (Recall, case B corresponds to a system 
where transitions occur on time scales that are large com- 
pared to the intrinsic time constant of the experimental 
probe (rp <^ {%}), while case D corresponds to a system 
where transitions occur on the time scale on the order 
of the intrinsic time constant of the experimental probe 
(rp « {%}).) The solution for both cases differ, which is 
solely due to difference in the magnitude of the rates of 
transition, see Tabled 

To garner insight into how the solutions for the auto- 
correlation functions depends on the magnitude of the 
rates of transition between biomolecular states, we now 
present two physically motivated guesses for the solution. 
We term these two guesses as the slow-switching limit 
and fast-switching limit solutions. (These solutions will 
be derived mathematically in Scct. lIVXI 'l If the molecule 
switches on time scales that is "long" compared to the 
intrinsic time constant of the experimental probe, then 
upon returning to is equilibrium position the probe will 
only observe a single biomolecular state. Upon averaging 
over many trajectories, switching will be infrequent and 
the probe should be capable of exploring each states sta- 
tistical properties. The auto-correlation function should 
then just be the average auto-correlation function of each 
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state 

{6xit)Sxm^{^e-'/^), (8) 
where explicitly 

i 

Here Pi is the probability that the molecule is in state i, 
ki is the associate spring constant of state i (Eq.[21) and 
Ti is the time constant of the probe when the biomolecule 
is in state i (ji = ^ /ki). This solution is termed the slow- 
switching limit solution, because the molecules switching 
(between states) is slow compared to the motion of the 
probe. On the other hand, if the molecule switches be- 
tween states on time scales that are very "short" com- 
pared to the intrinsic time constant of the experimental 
probe, then the molecule equilibrates between all states 
on time scales much shorter than those governing the 
motion of the probe. Therefore, the probe "sees" the 
average response of the molecule. The solution for the 
auto-correlation function should then correspond to that 
of an equilibrated system with the spring constant equal 
to the average spring constant of the biomolecule, 

{5x(t)5xm = (9) 

where 

(r-i) ^ {k)h (10) 

is the mean fc- weighted |40| inverse time constant. This 
solution is termed as the fast-switching limit solution, 
as the molecules switching (between states) is fast com- 
pared to the motion of the probe. In Sect. IIV Al we 
show that these solutions are indeed solutions for the 
auto-correlation function in the limit that the molecule 
switches either "slow" or "fast" compared to the motion 
of the probe. 

Figure |H1 displays these representations for the auto- 
correlation function. The slow-switching limit solution, 
Eq.|Sl is displayed as a blue line and the fast-switching 
limit solution, Eq.|51 is displayed as a red line. The auto- 
correlation function for case B lays on the physically mo- 
tivated slow-switching limit solution, Eq.|S| This is ex- 
pected, because case B corresponds to a system where the 
switching between the biomolecule states is "slow" com- 
pared to the decay time of the experimental probe. The 
time scales governing the switching for case D lies some- 
where in between the mentioned extremes. Therefore, it 
is reasonable to expect its solution to lie in between these 
two limits, as depicted in the figure. 

Although the solutions for the auto-correlation func- 
tion for both case B and case D differ, their separa- 
tion is only a few percent of their solution. From an 
experimental point of view it is important to under- 



stand the magnitude and origin of this separation, so 
as to determine whether or not an experiment can pro- 
vide sufficient spatial resolution to differentiate between 
various cases. To gain insight into how this separation 
may be influenced by the mechanical response of the 
biomolecule, we now expand our parameter space and 
consider two different cases, termed case B and case V. 
The parameters now governing the biomolecule (rates 
and spring constants) are the same as their counterparts 
(case B and D) except that we choose the spring con- 
stant for the biomolecule in state 1 to be larger, fc,„i = 
120 pN///m. (Recall, for case B and D kmi=2b pN//im, 
fcm2=5 pN/^m and A:p=20 pN//xm.) Although fairly 
larger, this new spring constant is still easily obtainable 
in muscle proteins [ll7lll0 llTlllll| . 

Figure 1^1 displays the results for the auto-correlation 
function in this new example. The separation between 
the two cases {B and D) is now much larger than the re- 
sults presented (B and D) in Fig. |S1 Additionally, there 
is a larger separation between the two physically moti- 
vated solution, represented by Eqs.|H|and|51 The increase 
is a direct consequence of the increase in the difference 
between the spring constants controlling the motion of 
the probe (difference between fci and k2, ki = k^^i + fcp). 
General features regarding the influence of the biomolec- 
ular properties on spatial separation is further explored 
in Sect.|V| 

Although there is an increase in the separation between 
various solutions the cousin of case B, case i3, also lies on 
the form expressed in Eq.|S| supporting the notion that 
this form is the proper solution for the auto-correlation 
function when the switching between biomolecular states 
is "slow" compared to the time scales governing the mo- 
tion of the experimental probe. As the rates increase 
(going from case B — > case V) the solution for the auto- 
correlation function approaches the solution represented 
by Eq.O similar to results found in Fig. |S1 This is con- 
sistent with the interpretation that the form of the auto- 
correlation function represented by Eq.|51 corresponds to 
the case where the biomolecule switches between states 
on times scales that are much faster than those governing 
the motion of the experimental probe. 

Although for these new cases the spatial resolution 
(difference between the two solutions) for the auto- 
correlation function has increased, limitations imposed 
by traditional methods (regarding time resolution) have 
not altered. Figure ^| displays the traditional averag- 
ing procedure presented earlier in this section. Again, 
the properties of the biomolecule when represented by 
case B can be adequately identified with this method, 
those represented by case V cannot. 

In this section we presented numerical solutions for 
the position of an experimental probe and its associated 
auto-correlation function, when the probe is attached to 
a biomolecule with fluctuating mechanical response. It 
was demonstrated that when the biomolecule fluctuates 
between states on time scales governing the motion of 
the experimental probe, sufficient temporal resolution 
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FIG. 9: Auto-correlation function for a system with a fluctu- 
ating spring constant. The parameters used to define the sys- 
tem here are similar to there counterparts displayed in Fig. |H| 
except the spring constant for the biomolecule in state 1 is 
120 pN//im. Case B and T) axe represented by blue circles 
and red circles, respectively. The solid lines refer to physi- 
cally motivated guesses, expressed in Eqs.|H| and |2l (blue line 
and red line, respectively) . 
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FIG. 10: Mean squared displacement of the probe {{x{t)^)T) 
averaged over some time window T, represented by gray dots. 
The top figure is the analysis for case B (T=30 ms) and bot- 
tom figure is for case T> (T=0.05 ms). Solid black line rep- 
resents the mean squared displacement of the probe in that 
particular state. 



is not provided for traditional analysis to be amiable 
for identifying the mechanical and kinetic properties of 
the biomolecule. Furthermore, we demonstrated that in 
this temporal regime the solution for the auto-corrclation 
function is affected by the biomolccular properties (rates 
and mechanical response) and can potentially be used in 



the analysis of experimental data. 

Our numerical study provided insight into the problem 
of interest. Although, they proved to be valuable in this 
context, numerical solutions for the equation of motion 
are not always the most efficient theoretical means to 
analysis experimental data; particularly when biomolcc- 
ular parameters of interest (rates and spring constants) 
are to be determined from the data itself (by least squares 
fit for example). For such cases analytic solutions pro- 
vide a more promising alternative. We now proceed with 
the mathematical formalism for calculating the auto- 
correlation function for a system fluctuating between me- 
chanically distinct states. Initially, we present bounds 
for the general solution to the auto-correlation function. 
Although, no general "closed" form solution for the auto- 
correlation function is known (at least to our knowledge) 
we present solutions in the form of two different series 
expansions, that are expressed around the two opposing 
bounds. Later, we present a general analysis regarding 
when the form of the auto-correlation function is appre- 
ciably influenced by the biomolccular properties. 



III. THEORETICAL FORMALISM: MOTION OF 
AN OVER-DAMPED PARTICLE IN A 
SWITCHING HARMONIC TRAP; BOUNDS ON 
THE AUTO-CORRELATION FUNCTION 

We now present our theoretical framework for calcu- 
lating the auto-corrclation function for a system that is 
in a thermal heat bath and is undergoing transitions 
between various mechanical states. In this section we 
prove that the two physically motivated solutions pre- 
sented in the previous section are indeed bounds to the 
general solution of the auto-correlation function. These 
bounds will determine the spatial resolution of the prob- 
lem at hand, by providing the possible range over which 
the auto-correlation function can fluctuate. In addition, 
these bounds serves are starting points for series solutions 
for the auto-correlation function presented in Sect. II VI 

To make this analysis more manageable (though not 
changing the principle of the calculation) we now assume 
that the equilibrium position is the same for all of the 
states and hence x°{t) = 0. In single molecule experi- 
ments this would correspond to measuring the motion of 
the experimental probe perpendicular to the orientation 
of the biomolecule, for example. The equation of motion, 
Eq.^ now becomes 

jx{t) = Fr{t) - k{t)x{t). (11) 

Additionally, we impose the stationary condition for the 
fluctuating spring constant. The stationary condition 
states that the average value of an arbitrary function 
of the fluctuating spring constant is independent of time 
translations. For a function G{k{t)) that is dependent on 
the spring constant at a single instant of time, the sta- 
tionary conditions states that its average is independent 
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of time, 

(G(fc(t)))=(G(fc)). (12) 

111 addition to the stationary condition we remind the 
reader that the assumptions presented in Sect. ^ still 
hold. These assumptions lead to the following impor- 
tant fact. Because the random thermal force is Gaussian 
and its first two moments (Eqs.|21and^ are independent 
of the state of the molecule, the random thermal force 
{Fr{t)) is independent of the state of the molecule and; 
therefore, statistically independent of the spring constant 
k{t). (This assumption does not hold if the damping con- 
stant was state dependent.) This results in a decoupling 
when averaging over functions of the random force and 
spring constant, 

(G(fc(0, • • • , k{t'))H{Fr{s), • • • , Fr{s'))) = (13) 
(GiHt),--- ,kit'))){HiF,is),--- ,F,is'))), 

where G and H are arbitrary functions. This fact will 
become important for the analysis of the auto-correlation 
function. 



A. Formal Solution for position of the probe 

To obtain bounds for the general solution to the auto- 
correlation function we first present a formal represen- 
tation for the solution to the position of the experi- 



mental probe. Equation 1111 is a linear, inhomogencous, 
stochastic, differential equation with stochastic coeffi- 
cients llSlSSISlElllllMEi- The general solution to 
Eg . 1111 (i.e. the position of the experimental probe) can 
be written formally as 

1 f°° ' 

x{t) = - / e-^o Ht-t")h dt"p^(^^ „ (^4) 
7 Jo 

where k{t) is the fluctuating spring constant. What lays 
difficulty into this problem is that the exponential is 
stochastic, due to the fluctuating spring constant. Evalu- 
ation of such a term and the associated averages arc quite 
diflicult, as we are aware of no g eneral solution. Only spe- 
cial cases have been solved |22[ 113, HE US • Although no 
general solution to this problem exists, the formal solu- 
tion for the position of the probe, Eq.^l can be used to 
obtain bounds on the auto-correlation function. 



B. Bounds on the auto-correlation function 

Using the formal solution for the position of the probe 
we now provide an upper and lower bound to the general 
solution for the auto-correlation function. These bounds 
correspond to the two physically motivated solutions we 
presented in the previous section, Eqs.|Hl and |^ Given 
the formal solution for the position of the probe. Eg . 1141 
the solution for the auto-correlation function is expressed 
as 



-1 i<oo noo I I 

((5a;(i)fa(0)) = I I dt'ds'e-fo '''"''^'-^''^^^e-'fo '''"''(-'")hFr{t-t')Fr{-s')j 

dt'ds'(|e-^o'''t"'=(t-*")/7e-/o'rf3"fc(-^")/7^^(i (15) 



T Jo Jo 

OO /'OC 



1 

7^ Jo Jo 



I 

We now exploited the fact that within our assumptions fluctuating spring constant. Eg . 11 31 The general solution 
that the random force is statistically independent of the to the auto-correlation function is now expressed as 



{Sxit)dx{0)) = — 



OO pOO 



T Jo Jo 

OO pOO 



T Jo Jo 



dt'ds'(^e- fo dt"k(t-t")he-SS' ''«"M-^")/7^2kBT7J(s' - {f - t)), (17) 



where the definition of the second moment for the random is 
force, Eq.^l has been used. Carrying out the integration 

over the dummy variable s' the auto-corrclation function ((5a-(t)(5x(0)) = (18) 



2kBT 

7 



OO 



dt'(^e-'^<>' dt"k(t-t")/'(^~ J^f-' ds"ki-s")h'' 
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We are now in position to obtain bounds for the auto- 
correlation function. To obtain an upper bound we start 
with the general definition of the auto-correlation func- 
tion, as expressed in and shift the last integral in 
the exponent by t' and change the dummy variable to t" . 
The general solution for the auto-correlation function is 
now expressed as 



{5x(t)5x{Q)) 



(19) 



2kHT r ^^// dt"Ht-t")h ^- dt"k{t' -t")h\ 



1 



To simplify our expression we define a function 
g{t,t',t"), which has the property 

75(M',t")^{^[;r^;;) . (20) 

The auto-correlation function may now be expressed with 
only one integral in the exponent, 



{6x{t)Sx{0)) 



2k„T 



7 Jt 



(21) 



For a given realization of the molecule over time (a par- 
ticular contribution to the average) an upper bound can 
be placed on the exponential 



-/„=^*'-'dt"g(t,t'.t") < 1 



2t' -t 



2t' -t 



^^1/ ^-{2t'-t)g(t.t' .t") ^ 

(22) 



This bound results because the arithmetic mean is 
greater than the geometric mean |3^ . This inequality 
leads to an upper bound on the auto-corrclation func- 
tion. 



{5x{t)5x{Q)) < 



(23) 



2k„T 



7 Jt 



dt'- 



2t' - t 7o 



2t'-t 



'dt"^e-(2*'-*)9(*^*'.*")^^ 



where we have taking the average inside of the inner in- 
tegral. To obtain our desired form for the upper bound 
we exploit the definition of g{t,t',t") (Fia. VM^i and the 
stationary condition (Ea. ll2|l . From the definition of 
g{t,t',t") the inner integral is expressed as 



r2t'-t 



dt"/p-(2t'-*)fe(t-t")/7 



r2t'~t 



dt" 



-{2t'-t)k{t'-t")/i 



(24) 



The stationary condition implies that the terms within following simplification 
the averages are independent of time, resulting in the 



2t'-t 



dt"(e-(2*'-*)M*-t")/7^ -t- / dt"(^e-(2*'-*)Mt'-t")/7^ 



rfi"/e-(2*'-t)fc/7\ 



2t'-t 



i'^g-(2t'-t)fe/7^ + {t' - t)(e-(2*'-*)'=/7^ 
(2i'-0(e-(2*'-*)"/^). 



^^„/ _(2t'-t)fc/A 



(25) 



Substituting this result into Eq.[531 the bound on the 
auto-corrclation function is now 

{Sx{t)Sx{0)) < ^ dt'(e-(2*'-*)fc/7y (26) 
Performing the integral we have 

{5x{t)SxiO)) < (^e-*/"^ (27) 



The right hand term is indeed the slow-switching limit 
solution for the auto-corrclation function. Proving that 
the slow-switching limit solution is an upper bound for 
the auto-correlation function. 



Having obtained an upper bound for the auto- 
correlation function we now proceed to obtain a lower 
bound. We now return to the general definition of the 
auto-correlation function, Eq.^l and express the fiuctu- 
ating spring constant in terms of its average value and 



13 



the fluctuations around that average, 

k{t) = (k) + Sk{t) (28) 

where 

6k{t) = k{t) - (k). 
I 



The average in Eq.lJHl is a sum over aU possible rc- 
aUzations of the molecule. For a given realization, the 
argument of the exponential is a fixed number, which 
allows us to exploit the inequality 

>1 + X, (30) 

I 



{Sxit)6x{0)) > ^^"dt'e<^)*/7e-2(fc)*77^1_ 



Exploiting the stationary condition, Eq.^J the later 
terms vanish 

{Sk{t)) = {k(t)) - (k) 
= 0. 

Therefore, the lower bound expressed in Eq.|^ can be 
written as 

{6x{t)Sx{0)) > ^ /"dt'eW'/^e-^W'VT. (32) 
T J t 

Performing the integrals we get our desired result 

{Sx{t)Sx{0)) > ^e-*<""'>. (33) 

The lower bound is the fast-switching limit solution and 
has the physical interpretation that due to the fluctua- 
tions in the spring constant over time, the displacement 
of the probe is more correlated to past events than an 
equivalent system with a single spring constant equal to 
the average spring constant of the fluctuating system. 
The auto-correlation function not only reflects the av- 



Due to the stationary condition the average value of the 
spring constant is time independent. With this new ex- 
pression for the fluctuating spring constant, the auto- 
correlation function is now expressed as 



(29) 

I 

for arbitrary X. With this inequality a lower bound can 
now be placed on the auto-correlation function 



dt"5k(t-t")h + j'^ ' ds"5k{~s")hyi 

[j\t"(dk{t-t"))h+ 1^ ' d.s"(dk{-s"))/j)y (31) 

I 

erage response of the biomolecule but also correlations 
between the biomolecular states. 

The fact that these two solutions bound the general 
solution to the auto-correlation function, 

^^e-*<-^> < {Sx{t)Sx{0)) < {^e-'/^), (34) 

also appeals to physical intuition. Consider the fanci- 
ful yet instructive case where the magnitude of the rates 
for the biomolecule can be tuned, similar to the analysis 
presented in Sect.^ Consider a parameter a which de- 
scribes this tuning, {kij} — > {aKij}, for i ^ j. If we ini- 
tially set a <C 1 then the solution for the auto-correlation 
function will correspond to the slow-switching limit. As 
the parameter a increases, switching occurs more quickly 
and the solution for the auto-correlation function devi- 
ates from the slow-switching limit solution and begins to 
approach the solution in the fast-switching limit. As the 
parameter further increases, such that a ;3> 1, we arc in 
the fast-switching limit. (This type of tuning is displayed 
in Fig. Inland 1^) As we tune this parameter a the solu- 
tion for the auto-correlation function starts at one limit 
and approaches the other, always lying in between these 
two bounds. 



{Sx{t)Sx{0)) 



2k^T 

1 



dt'e^''>*/'^e-2W*'/7/g-(Jo' dt"'5fc(t-t")/7+/o'"' 'is"'5fe(-s")/7) 
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In this section we have obtained an upper and lower 
bound for the general solution to the auto-correlation 
function. These bounds were found by exploiting prop- 
erties of the formal solution, However, the value 
of the formal solution as a tool amiable to the analysis of 
experimental data is less evident. Transforming the for- 
mal solution into a tractable "closed" form proves to be 
difficult, as we arc aware of no such solutions. Solutions 
more tractable for the analysis of experimental data can 
be sought by a series expansion |^|33,0- ^^^^ sense, 
we follow earlier developments, though our method dif- 
fers from previous methods in that the series is expanded 
around two important experimental limits: the fast- and 
slow-switching limits. 

A great benefit of presenting these two series solutions 
is that they approach the exact solution from oppos- 
ing bounds. To generate these two series solutions we 
choose to first present a single series solution for the auto- 
correlation function, upon which both limiting forms can 
be derived. Upon deriving the limiting forms we trans- 
form the single series into the two desired series solu- 
tions. Generating the two series solutions in this fashion 
will lead to a robust analysis of the error associated with 
truncating the series at a given order (Sect.EJ. Giving 
a prediction regarding the convergence of each series a 
priori. 

IV. SERIES SOLUTION FOR THE 
AUTO-CORRELATION FUNCTION 



We now present a series solution for auto-correlation 
function, upon which the two series solutions which ini- 
tiate at the limiting forms will be derived. To determine 
this first series solution we express the solution for the 
position of the probe x{t) itself in the form of a series. 
Given the solution for the position of the probe x{t), 
the auto-correlation function can be determined. The 
position of the probe x{t) satisfies the equation of mo- 
tion represented in Eg . 11 II and obtaining its solution is a 
formidable task as the equation of motion depends on two 
stochastic variables: Fr{t) and k{t). To proceed we find 
it convenient to reformulate the equation of motion into 
a slightly different form, whose utility will be revealed 
below, 

-fX=: Fr{t) -kaX + f{t)kdX. (35) 

Now the switching between different mechanical states is 
captured in the presence of the function f{t). By writing 
the dynamics in this way, we are required to introduce 
new notation. In particular, ka is the equally- weighted, 
average spring constant, averaged between the maximum 
and minimum values, 

— (^max ^ ^min)/2, 

kd is the equally-weighted, difference spring constant. 



which reflects the difference between the maximum and 
minimum values, 

— (^max ^min)/2, 

and f{t) is a random variable which determines the spring 
constant of the system 

k ~ k- 

fit) = -—, if in state i. (36) 

kd 

We note that the choices for ka and kd are arbitrary 
though they have the advantage that this choice bounds 
the magnitude of the fluctuating variable f{t) (|/(t)| < 
1) and the ratio of the spring constants {\f{t)kd/ka\ < 1) 
which provides analytical advantages. 

The new formulation can be interpreted as follows. 
The probe is in a harmonic trap with spring constant 
ka, subjected to two random forces: Fr{t) mimics the 
thermal heat bath and f{t)kd emulates the switching of 
the states of the molecule. The form of the switching 
term is now a small perturbation when compared to the 
harmonic trap ka, \fit)kd/ka\ < 1. This switching term 
"kicks" the probe either toward the origin or away from 
it, depending on whether state i has a spring constant ki 
greater than or less than ka, respectively. 

To obtain a scries solution for the auto-correlation 
function we first promote the fluctuating spring term, 
-~f{t)kdx{t), in Eq.|3Slto the inhomogeneous part of the 
differential equation (treat as an external force). Next 
standard solutions to the problem can be used to for- 
mally rewrite the solution to Eg. 1351 as 

x{t) ^ Xrit) +Xf{t). (37) 

The first term Xr (t) is the displacement induced by ther- 
mal fluctuations in a harmonic well described by spring 
constant ka, namely, 

1 Z""" t' 
x^{t) = - dt'e^^Frit-t'), (38) 

7 Jo 

where Ta ~ l/ka is the equally weighted average time 
constant. The latter term Xf{t) is the displacement in- 
duced by the fluctuating spring f{t)kd in the same har- 
monic well described by spring constant ka, namely, 

1 ' 
Xf{t)^-l dt'e-^ f{t~t')kdx{t-t'). (39) 
7 Jo 

The motion of the probe is now described by two types of 
random "kicks" . The thermal heat bath induces random 
"kicks" described by motion Xr{t) and the switching of 
the macromolecule induces random "kicks" described by 
motion Xf{t). 

The solution for x{t) is now a functional of itself (due 
to the X f (t) term) and can be solved itcratively in terms 
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of the displacement due to thermal motion: 

x{t) = (40) 

Y.\ — ] / n '^^^ e-^FrXt, {Tj})Xr{t - T„), 

where = j/kd is the equally weighted difference time 
constant and 



r„ = h + t2 + h 



t,. 



(41) 



^„(t,{T,}) = /(f-Ti)/(f-T2)---/(t-r„);(42) 
and the following conventions are used: 

^o(i,{T;})^l 
and for arbitrary function i/(t, {Tj}) 

poo 

/ l[dUH{t,{T,})^H{t). 

Following these conventions the zeroth order term in 
I 



Eq.0U|is simply Xr{t), which is just the thermal motion. 
Addition of the next term provides an approximation for 
the solution, whose form is 

1 1"°° I 

x{t)KXr{t) + - / dt'e-^ /^'^f{t-t')kdXr{t-t'). (43) 
7 Jo 

This solution has the physical interpretation that the mo- 
tion of the probe is a linear combination of the thermal 
motion Xr{t) plus the thermal motion subjected to the 
fluctuating spring force ~f{t)kdXr(t). The higher or- 
der corrections describe correlations which impose self- 
consistency for x{t). Truncation at second order would be 
reasonable if thermal motion were dominate over switch- 
ing motion, Xr{t) ^ Xf{t). In general, however, we 
are not interested in such cases, because extraction of 
the molecular properties from experiments becomes quite 
difficult when thermal noise dominates over switching. 
Moreover, we find it more useful to expand the auto- 
correlation function around two physically important, 
limiting cases. 

Using Eq.0U|a series solution for the auto-correlation 
function is obtained: 



{Sx{t)6x{Q)) = {-) / n n dUdU,e-^^{T„{0, {T,})T,-,,{t, {Tj,})xr{'T„)xr{t - T^')). (44) 



Using the fact that within our approximations the fluctu- 
ating spring constant is statistically independent of the 
thermal force, Eq.^l dictates that the motion of the 
probe induced by thermal fluctuations Xr{t) is indepen- 
dent of the switching function f{t). The averages over the 
thermal induced displacements Xr(t) and the stochastic 
variable f{t) then decouple, 

{fitl)---f{t^)Xr{t)Xr{t')) - {fih)---f{t„,)){xr{t)Xrit')) 

= ^e-^(/(ii)---/(U). 



The auto-correlation function may now be written as 



ri.Ti'— t—li'—l 



\i + Tn-T ,\ + Tn + T 



(45) 



r 



where the stochastic features regarding the switch- 
ing of the molecule is embedded in the term 
{J^n{0,{Tj})^n'{t,{Tj'})}- Although Eq.EHl provides 
a valid series for the general solution for the auto- 
correlation function, we find it more illuminating and 



beneficial to begin such a series from the two obtained 
bounds found for the auto-correlation function. We now 
derive (mathematically, as opposed to physically) the so- 
lutions for the auto-correlation function in the fast- and 
slow-switching limits directly from Eq.021 This deriva- 
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tion will allows us to express the general solution for the 
auto-correlation function in the terms of two series solu- 
tions initiated at these limiting forms. 



A. The fast- and slow-switching limits 

In our previous analysis, Sect. Ill 51 we introduced two 
possible solutions to the auto-correlation function based 
on physical arguments and proved in Sect. Illll that these 
two solutions provided an upper and lower bound to the 
general solution for the auto-correlation function. We 
denoted these two limits as the fast-switching limit and 
the slow-switching limit. In this section we derive these 
solutions starting from the series solution expressed in 
Eq.EHl 

Mathematically, the fast-switching limit is defined as 
the case when Kg 3> 1/tq, where Kg is the slowest rate 
of transition between the different states (ks = min(Kij) 
for i ^ j). This states that the molecule switches be- 
tween states on times scales that are much shorter than 
those governing the motion of the probe. The slow- 
switching limit is defined as the case when K{ <C I/tq, 
where Kf is the fastest rate of transition between dif- 
ferent states {k{ = max(Ky) for i ^ j). This states 
that the time scales governing the motion of the probe 
is much shorter than those governing the transitions of 
the molecule. For the problems of interest Tq (where 
Ta = 7/(fcp + (fcr'r" + C'")/2) and - niax(fc„„) 

and fcJJJ™ = min(fcnii)) can be considered as a measure 
of the correlation time of the experimental probe, as the 
compliance of the probe (1/fcp) is generally much smaller 
or of order of the compliance of the biomolccule (1/fcmj) 
in any state. (Sec discussion in Appendix IbI) 

The auto-correlation function, Eq.^^l contains aver- 
ages over the random variable f{t) along with a decaying 
exponential with time constant Ta- In the fast-switching 
limit the averages over the random variable f{t) decouple 
on time scales much faster than any other in the problem. 
Therefore, we can assume that any correlation between 
the random variables (/(t)) at different times are negli- 
gible and make the following approximation: 

(^„(0,{T,}).F„.(i,{r/}))= (46) 

(/(-Ti) • • • /(-r„)./(t - TiO ---fit- « (./)"^"', 

where 




Here fi is the value of / when in state i, Eq.lJHl and we 
have used the stationary condition for f{t), Ea.lT^ 

In the slow-switching limit the averages over the ran- 
dom variable f(t) will remain nearly constant on time 
scales where the integrand is non-negligible. Therefore, 
we can assume that the averages have no dependence on 



time differences and we make the following approxima- 
tion 

{j^„{OAT,Wn'{t,{Tr}))= (48) 

(fi-Ti) ■ ■ ■ f{-T,,)f{t - Ty) ■■■fit- r„o) « (/"+"'), 

where 

i 

An important result is the fact that in these limits the 
averages over the random variable are independent of any 
time variables. It is shown in Appendix lO that if 

{f{tl)f{t2)...f{tm)) - C„,, (50) 

where Cm is independent of the time variables, then the 
integrals in the series fEa. l45|l can be done explicitly and 
the solution for the auto-correlation function may be ex- 
pressed as 

,rp oo / \m rn / 4- \ j 

(MWO)) = ^e-/.^ ^ (I) a.E^(f) ■ 

(51) 

For the two limiting cases (Eqs.^l and I49|) this result 
can be summed in closed form and the resulting auto- 
correlation functions arc: 

r ^e-*<--^> Kg » IK, 
{5x{t)5x{Q)) = I (52) 
i (i^e-*/-) «f«lK. 

These are precisely the physically motivated solutions we 
have presented in the previous sections. 

As expressed in Sect. IIIBI the rate of switching influ- 
ences the effective interaction between the biomolccule 
and experimental probe. In the fast-switching limit the 
biomolccule switches between states so "fast" that it 
equilibrates on the time scales shorter than those gov- 
erning the motion of the probe (conveyed mathemat- 
ically in Ea. l46|) and; therefore, the probe "sees" the 
biomolccule as a stable single-state system with an effec- 
tive spring constant equal to the average spring constant 
of the biomolccule (km), 

(A:) = EP,fc, = A:p + (fc„,). (53) 

i 

In this case the motion of the probe cannot be used to 
delineate the different states, but only infer an effective 
interaction. In the slow-switching limit thermal decay 
causes the experimental probe to loose memory of any 
previous position on times scales that are much smaller 
than the lifetime of each state (conveyed mathematically 
in Ea. l48|) and; therefore, the probe is able to fully resolve 
each states auto-correlation function, properly weighted. 
It is reassuring to us that these limiting cases appeal to 
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physical intuition: in the fast-switching hmit it is the av- 
erage spring constant that is physically observed, while in 
the slow-switching limit it is the average auto-correlation 
function that is physically observed. 

B. Series expansion around limiting forms 

Armed with the derivation of the two limiting forms 
from our general scries expression, we now express the 
general solution for the auto-correlation function in the 

I 



^ ' n,n'—0 ^ ^ 



where 

^'Sl^.mi.m-'}) ^ {Tr.io,{T,})T,At,{Tr}))-{fy'^ 

and the latter term ((/)"^" ) is the contribution to the 
I 



^ ' n,n'— 



where 

-Dl'U^^mATf}) ^ (.F„(0,{r,}).F„-(t,{T,,}))-(/"+«') 

and the latter term ((/"^" )) is the contribution to the 
limiting form. Here the superscript (0) denotes that the 
expansion is around the slow-switching limit {nfTa <C 1). 
Note, that by definition 

^0,0 ~ ^0,0 — ^0,1 ~ -^1,0 ^ ^0,1 ~ -^1.0 — '■^ y^^J 

and the leading order correction is second order in both 
series. 

The two series, Eqs.|^ and |55l are general solutions 
to the auto-correlation function which initiate at oppos- 
ing limits and approach the solution from opposite direc- 
tions, adding to utility of expanding the solution around 
these two limiting forms. The expansion terms in Eg . 1541 
must add up to a positive number and those in Eg . 1551 
must add up to a negative number. 

When the solution for the auto-correlation function 
corresponds to one of the limiting forms the auto- 



form a series initiated from either of the limiting forms. 
The deviation of the auto-correlation function from the 
limiting forms reveals the time scales on which the 
biomolecule fluctuates between various states, as was il- 
lustrated in Sect . ITlBl fFig. 151 and 1^ . To obtain the two 
new general series solutions we now add to the general 
solution, Eq.^ni the limiting forms, Eq.E3 and subtract 
from each term the contribution that is already captured 
in that limiting form (Eg .1471 or I49|l . 

Initiated at the fast-switching limit the general solu- 
tion for the auto-correlation function takes the form 



J] n '^^^^^^■^i «' {^^}' iTr})e ^ , (54) 

i=i i>=i 

I 

limiting form. Here the superscript (oo) denotes that the 
expansion is around the fast-switching limit (ksTq ^ 1). 
Initiated at the slow-switching limit the general solution 
for the auto-correlation function takes the form 



n fl dUdtMun'it, {Tj}, {T,v})e-"^"""""'r""^""', (55) 

i=l i' = l 

I 

correlation function cannot be used to convey the rates 
of transition, as these limiting solutions only depend on 
probabilities. The magnitude of the rates of transitions 
are conveyed through the terms T>1^, and X)^"^, . When 
the time dependence for these terms reveal themselves 
through the solution to the auto-correlation function can 
the magnitude of the rates of transition be determined. 
This time dependence will be revealed if these terms 
(2?'°°-! and p'"',) vary on the time scale t,,- We there- 
fore, expect that measurements of the auto-correlation 
function would be most applicable to determine the rates 
of transitions when the rates are of order of the inverse 
time constant of the experimental probe, {Ky} ~ l/^a, 
for i ^ j (not near either of the limiting forms as exem- 
plified in Fig. IHlandini). This is a time resolution that is 
inaccessible by traditional analysis such as that presented 
in Sect.ini 

Although we believe that this method will be most ami- 
able to experimental analysis when the solution for the 
auto-correlation function is not near one of the limiting 
forms, analysis of the auto-correlation function may still 
provide valuable information regarding the rates of tran- 
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sition when this condition docs not hold. For the case 
when the system is in the fast-switching hmit, analysis 
of the auto-corrclation function does not provide enough 
temporal resolution to resolve the absolute rates of tran- 
sition, the probe "sees" the average spring constant; how- 
ever, if the system is in the fast-switching limit, analysis 
of the auto-correlation function provides a lower bound 
on the rates, as by definition {nij} ^ 1/Ta for i ^ j. 
This lower bound is far greater than that which can be 
concluded by traditional methods, as that presented in 
Sect.|lll(see Fig.H(c) and (d)). When the solution for the 
auto-correlation function is properly represented by the 
slow-switching limit traditional methods like that pre- 
sented in Sect. |n] are applicable and far superior. 

If the calculation of the auto-correlation function from 
the series expansion (Eq.|^ or Ea. l55|l and its compar- 
ison with experimental results is practical, the synopsis 
presented here supports the notion that the calculation 
of the auto-correlation function may be used as method 
to determine the rates of transitions for macromolecular 
systems, those which are inaccessible with current ana- 
lytic methods. To help garner insight into the applica- 
bility of this method we now explore issues regarding the 
convergence of both scries to the general solution. 

V. CONVERGENCE AND APPLICABILITY 

A key to the applicability of this method will be on 
the ease at which the series expansions converges to the 
general solution; in this section we present an analysis 
regarding this convergence. We then offer a concrete un- 
derstanding of how the biomolccular properties affects 
both the convergence of each series onto the general so- 
lution and the applicability of this method to study rates 
of transition by focusing on the special case of the two 
state problem. The two state problem is further explored, 
expressions for each series expansion (up to third order) 



are presented for a Markov process, expressions which are 
applicable to experimental analysis. Guided by the in- 
sight obtained in the section, regarding the convergence 
of each series to the general solution, we then test these 
expressions to predict the form of the auto-correlation 
functions calculated numerically in Sect. ITlBl 

We now generate expressions which lend insight into 
the convergence of each series. We portray the conver- 
gence of both series through the analysis of a new series, 
which we term as the difference series. This difference 
series measures the separation between both series so- 
lution when truncated at a given order. The utility of 
analyzing the difference series is that it can be written in 
a simplified closed form and; therefore, readily studied. 

Although this difference series does not provide the 
error associated with truncating a particular series, it 
stills proves to be a practical measure for the conver- 
gence of each series onto the general solution. This is 
because both series start at opposing bounds and must 
converge to the same result. In addition the utility of 
this methodology for determining the rates of transition 
will depend on the condition that the solution for the 
auto-corrclation function is separated from both limit- 
ing forms and; therefore, "equally" separated from both 
series. The separation between both series at a given 
order should then serve as a reasonable upper-bound es- 
timate |4lj for the separation of an individual series from 
the true solution. 

We define the difference series as 

AA^it) ^ {5x{t)5xi0))^^^ - {Sxit)SxiO))l^\ (57) 

where A.Am (t) is the difference between both scries when 
they are truncated at order m. The superscript on the 
correlation functions denotes the particular series from 
which the auto-correlation function is inferred and the 
subscript denotes the order of truncation. This differ- 
ence, at the m'^'-order, is 



^ ' \ / n,n'—0 ^ ^ ^ 2—1 I' — l 



where the sum is restricted to n + n' < m. Noting that 

((/)"^" — ^/"+" ^) is independent of time wc make use 
of Eg . 1511 and rewrite this difference as 



keT 



1 



't/r 



n=0 ^ 



1 



1 



(58) 



((/)" - (/")) E 



J! 



The difference series is now expressed in a simplified 
form, more amiable to analysis. The difference series 
can be further simplified for the t = Q case, where the 
difference series can be summed into a closed form 



AA40) 
k„T 



k/ ka 



[l-{k)/ka 
(k) 



,ni+l 



(59) 

The advantage of exploring the convergence of the differ- 
ence series AAm{t)j as oppose to each individual series, 
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is that it is much simpler to evaluate. The difference se- 
ries only depends on equilibrium averages ((/™)) instead 

of integrals containing correlations (2?^°^' , T>1^\, ) . 

For large times the difference series will be dominated 
by the decaying exponentials, even though the latter 
summation [30j (j-summation in Eg . 158(1 increases with 
time. Therefore, the magnitude of the difference series 
will be dominated by short times, t < Ta- Because of this 
we use the simple < = difference, Eq.EHl as a measure 
of the error in either series; which is more amiable to 
analysis. 

As an aside, we note that if expansion terms are cal- 
culated for one series expansion then the corresponding 
terms in the other series expansion can simply be calcu- 
lated through the employment of Eqs.EI a-nd [SSI This 
can greatly reduce the work in calculating the second 
series. 

Using the difference series as our guide, we now dis- 
cuss how the convergence of the series depends on var- 
ious parameters in the problem and determine a range 
for which this method may be applicable for the analy- 
sis of experimental data. For concreteness we focus on 
the two state problem, where biomolecules exist in two 
distinct mechanical states. We envision an experimen- 
tal situation similar in fashion to that in Fig. ^ and [3| 
(except where the macromolecule can only access two 
states). The two state problem is both revealing and 
serves as the basis of a number of actual applications of 
these ideas. In addition to determining a range of param- 
eters for which this method may be applicable we present 
truncated series solutions (up to third order) for the two 
state Markov process; giving functional form to the so- 
lution to the auto-correlation function when it deviates 
from the limiting solutions. 



A. The two state problem 

One of the key ambitions of this section is to determine 
a range of biomolecular parameters for which the analy- 
sis of the auto-correlation function will be beneficial for 
revealing such properties, particularly the rates of transi- 
tion. As noted, the form of the auto-correlation function 
depends on the rates of transition more sensitively when 
the rates are of order of the inverse of the time constant of 
the experimental probe, {Ky} ^ 1/Ta,i 7^ ]■ Under such 
conditions the form of the auto-correlation function is in- 
fluenced by the magnitude of the rates of transition and 
not just the probabilities of lying in a particular state. 
Therefore, a criterion for applicability is that the rates 
of transition are of order of the inverse time constant of 
the experimental probe. 

However, satisfying this criterion alone will not guar- 
antee the applicability of this method. For an experimen- 
tal study to exploit this method it must be able to at- 
tain sufficient spatial resolution to differentiate between 
various possible solutions, that is it must be capable of 
discerning between solutions in the fast-switching limit. 



the slow-switching limit or somewhere in between. The 
spatial resolution is directly correlated with the separa- 
tion of the two bounding solutions, the fast- and slow- 
switching limits; a smaller separation requires greater 
spatial capabilities of the experimental apparatus. As we 
have demonstrated in Sect. irTBl the separation between 
the slow-switching and fast-switching limits is affected by 
the ratio of the spring constants between the two states, 
as exemplified in Fig. [S] and |51 Imposing applicability 
will impose a restriction regarding the ratio of the spring 
constants of the states. Additionally, we will show that 
applicability will enforce limits on the probability of the 
system to lie within a given state. This later enforce- 
ment is simple understood from the fact that if one state 
dominates, the influence of the other will not apprecia- 
bly affect the form of the auto-correlation function and; 
therefore, difficult to detect. Finally, for this analysis 
to determine the rates of transition from experimental 
measurements of the auto-correlation function requires 
terms within the series to be computed. Calculations of 
such terms can be tedious, leading to the final constraint 
that applicability requires adequate convergence to occur 
within a limited number of expansion terms. This last 
condition combined with the need for deviation from the 
limiting forms restricts the parameter space for which 
this method should be applicable. 

We have given a number of conditions, which will 
restrict the parameter space, for when we believe this 
method will be applicable to experimental analysis. To 
determine this space of molecular parameters we now 
proceed by going through a number of refinement proce- 
dures. The first two conditions we explore are concerned 
with the requirement that to amiable to experimental 
analysis (for determining the rates of transition) requires 
the solution for the auto-correlation function to be "well" 
separated from both limiting forms. This will entail a 
condition regarding the probability of the system to lie 
within a given state and a condition regarding the ratio 
of the spring constants. Finally, in order to limit the 
number of integral evaluations (number of terms in the 
scries) required for a desired convergence, an additional 
condition regarding the spring constant ratio is imposed. 

1. Parameter space of applicability for the two state 
problem 

As mentioned at the end of Sect. IIVBI the auto- 
correlation function can in principle be used to deter- 
mine the individual rates of the macromolecule when 
it differs from both limiting forms (((5a;(t)(5x(0))p°°'' or 

((5a:(i)(5a:(0))^°''). Because both limiting forms provide 
opposing bounds to the solution for the auto-correlation 
function, a necessary condition for applicability requires 
the bounds to differs from each other more than the avail- 
able spatial capabilities of the experimental apparatus. 
Therefore, to obtain an initial range of parameters for 
which this method may be of use, we study the differ- 
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ence series at t = and m = (A^o(O)), which serves as 
a measure of "nearness" for both limiting forms. (AAo(O) 
is the difference between the hmiting forms.) When 
AAo(O) is large then the limiting forms are well sepa- 
rated and the solution for the auto-correlation function 
is not "near" the limiting forms. (Recall, we are consid- 
ering cases for which {Ky} « I/tq,^ ^ j and; therefore, 
the solutions lies in between the limiting forms.) 

Prior to proceeding with this analysis, we need a defi- 
nition of "nearness" . For concrctcness we consider an ex- 
periment whose error tolerance for calculating the auto- 
correlation function is of the order of 5% of the mean- 
squared displacement of the experimental probe (a;^). To 
obtain a range of parameters we seek situations where 
both limiting forms differ by more than 5%. 

Figure ^] displays the difference series at i = and 
m = 0, normalized by the solution in fast-switching 
limit. (The fast-switching limit is an approximate mea- 
sure of (x^).) The results are plotted as a function of 
the spring constant ratio, at various probabilities. The 
limiting forms give similar answers (A^o(O) is small) if 
one state is dominate (very large probability), so long 
as the spring constants do not vary by many orders of 
magnitude. Here, we only consider experimental situa- 
tions where the spring constants in both states are of the 
same order of magnitude. Cases for which this differs are 
unlikely for the experimental situations of interest. 
Therefore, to consider cases where the difference in the 
bounds is greater than 5%, a necessary condition that we 
impose is the probability of being in either state must be 
greater than 10%. 

Having produced an initial bound on the probabilities 
for applicability of this analysis, we now seek an initial 
lower bound on the ratio of the spring constants. For 
argument sake we now choose state 1 to be defined such 
that fci > k2. When fci w fc2, either the mechanical 
properties of the macromolecular states are similar or the 
spring constant of the experimental probe is dominate 
over that of the molecule. In such cases the t — Q solu- 
tions for both bounds are nearly identical and it should 
become quite difficult to use the general solution for the 
auto-correlation function to extract the rates of transi- 
tion. For the bounds to have a difference greater than 
5% (our predefined tolerance) then the spring constant 
ratio (fci/fc2) should be greater than 1.5. This puts a 
limitation on the stiffness of the experimental probe as 
compared to that of the macromoleculc. For the spring 
constant ratio to be greater than 1.5, then the macro- 
molecular stiffness in state 1 must be greater than 1/2 
the stiffness of the experimental probe. 

Satisfying these two conditions, the solution for the 
auto-correlation function deviates from the limiting 
forms on the order of 5% or greater. However, this does 
not guarantee applicability. The next condition that we 
impose regards limiting the number of terms in the series 
needed for adequately convergence to the desired solu- 
tion. Here we are interested in cases for which the series 
solutions, Eas.l54l and 1551 are within the predefined tol- 




FIG. 11: Difference series (t = 0, m = 0) Aj4o(0) normalized 
by the solution in the fast-switching limit {x^)ao as a function 
of the spring constant ratio, at various probabilities. (Pi is 
the probability of being in state 1.) Note that when this ratio 
is 1 both statse are mechanically indistinguishable. 



erance (5% of {x'^)) of the exact solution only after the 
calculation of three expansion terms. Figure iniplots the 
t = difference series AAm{Q) (normalized by k^T /ka) 
for two different spring constant ratios at various proba- 
bilities, as a function of the number of expansions terms. 
We choose the spring constant ratio to be 2 in figure 
(a) and 3 in figure (b), recalling that it must be greater 
than 1.5 to be within our initial refined parameter list. 
From Fig. I12f a'l it is evident that truncation of the series 
solution at second order is sufficient to allow for conver- 
gence to within 5% accuracy of the exact solution. From 
Fig. I12r b'). however, it takes more terms. To be within 
the experimental tolerance of 5%, the series must contain 
at least four terms. 

As stated in Sect. |llll (Eq.]^ there are no first order 
corrections; therefore, only one (three) expansion term 
(terms) are needed to be within this predefined error tol- 
erance for these cases. For spring constants ratios greater 
than 3 more terms are needed. Therefore, this analysis 
gives an upper bound on our spring constant ratio to be 
less than 3. 

The analysis presented in this section gives a range of 
parameters for which we believe this theory should be 
useful and practical for determining the rates of transi- 
tion (for the two state problem) from experimental mea- 
surements, those which are not accessible from tradi- 
tional means. To determine the rates of transition the 
solution for the auto-correlation function must deviate 
from both limiting forms, requiring the rate of transitions 
to be of the order of the inverse time constant of the ex- 
perimental probe, {Kij} « l/Ta,i ^ j. To be within 5% 
accuracy of the exact solution and not near one of the lim- 
iting forms, the probability of being in either state should 
be greater than 10% and the ratio of the spring constants 
should also be greater than 1.5. To limit the number of 
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2. Solutions for the Markov process 
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FIG. 12; Difference series Am{0) normalized by (A;)/kBT 
plotted as a function of the number of expansion terms, for 
k\/k2 = 2 (a) and k\/k2 = 3 (b), for various probabilities. Pi 
is the probability of lying in state 1. 



Our system of interest is a two state Markov process, a 
system consisting of tvi^o distinct mechanical states (dif- 
ferent spring constants) whose transitions are local in 
time (independent of all previous configurations). Ap- 
pendix ^ presents the detailed mathematics describ- 
ing the two state Markov process and the solution for 
the auto-correlation function truncated at third order, 
Eas. lDQI and lDllI Here we compare those solutions with 
the numerical solutions presented in Sect. Ill Bl and deter- 
mine the degree of convergence. 

Figure El redisplays the numerical solutions for case D 
presented in Sect. IIIBI Recall, case D corresponds to a 
system where the rates of transition are of order of the 
inverse of the intrinsic time constant of the experimental 
probe (see Tabled). In addition, solutions correspond- 
ing to the fast-switching limit (solid blue line) and slow- 
switching limit (solid red line) . see Fig. |S1 are also redis- 
played. To make connections with the previous analysis 
regarding convergence, the ratio of the spring constants 
for the two states is 1.8 (fci = 45pN/^m, fc2 = 25pN/^m) 
and the probability of lying in state 1 is 30%. These 
conditions nearly satisfy those displayed in Fig. ^1 (a) 
and close to the limit we imposed for applicability. From 
Fig. El (a) we expect the convergence of both series to be 
within 1% of the exact solution. The analytic solutions 
truncated at third order are also displayed in the fig- 
ure. The series expressed around the fast-switching limit, 
Eg. lDQI is represented by the red dashed-line and the se- 
ries expressed around the slow-switching limit, Ea. lDllI 
is represented by the blue dashed-line. As expected, both 
series truncated at third order accurately portrays the 
numerical solution for case D. 



calculated terms in the series to three or less, the ratio 
of the spring constants should also be less than 3. 

This range of parameters is somewhat limited; how- 
ever, this method still provides an alternative to stan- 
dard analysis, expanding the temporal capabilities al- 
loted in experimental studies. Having obtained this range 
of parameters for applicability we now present calculated 
terms in both series, up to third order, for the two state 
Markov process. These solutions are amiable to exper- 
imental analysis, whose form depends on the value of 
the individual rates. Guided with our understanding of 
convergence obtained in this section, we then compare 
our truncated analytic series solutions to the numerical 
results presented in Sect. IIIBl 



The molecular parameters expressed in case D are near 
the region of minimal acceptable spatial resolution. Spa- 
tial resolution was enhanced when the spring constant as- 
sociated with state 1 was increased to ki = 140pN/^m, 
which we termed as case T>. This results in a spring 
constant ratio of fci/fc2 = 5.6, which is somewhat out- 
side the predicted bounds for applicability. To test this 
idea, Fig. El presents the numerical solutions for case V 
along with the truncated series solutions. The numer- 
ical solutions are again represented by red circles, the 
solutions in the fast-switching limit is represented by the 
solid blue line and slow-switching limit is represented by 
the solid red line. The analytic solution expressed around 
the fast-switching limit, Ea. lDQI is represented by the red 
dashed-line and the analytic solution expressed around 
the slow-switching limit, Ea. lDllI is represented by the 
blue dashed-line. As expected, both series truncated at 
third order do not accurately portrays the numerical so- 
lution for case T). 
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FIG. 13: Auto-correlation lunction for a system with a fluctu- 
ating spring constant. Numerical solutions for case D is rep- 
resented by the red circles. The solid red line is the solution in 
the fast-switching limit and the blue line is the solution in the 
slow switching limit. Series solutions for the auto-correlation 
function, truncated at third order, are represented by the red 
dashed-line when initiated from the fast-switching limit solu- 
tion (Eg. lDQt and by the blue dashed-line when initiated from 
the slow-switching limit solution (Eg . ID 111 . 
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FIG. 14: Auto-correlation function for a system with a fluctu- 
ating spring constant. Numerical solutions for case T> is rep- 
resented by the red circles. The solid red line is the solution in 
the fast-switching limit and the blue line is the solution in the 
slow switching limit. Series solutions for the auto-correlation 
function, truncated at third order, are represented by the red 
dashed-line when initiated from the fast-switching limit solu- 
tion (Eg . ID9t and by the blue dashed-line when initiated from 
the slow-switching limit solution (Eg . ID 111 . 



VI. CONCLUSION 

111 this work we presented an analytic method applica- 
ble to the analysis of data from single-molecule switch- 
ing experiments. This method regards the analysis of 
the auto-correlation. The auto-correlation function was 
taken as the vehicle for analysis, because it varies on 
time scales governed by the intrinsic time constant of the 
experimental probe and; therefore, its functional form 
should be effected by biomolecular transitions which oc- 
cur on such corresponding time scales. This is an im- 
portant attribute, because we presented arguments (the- 
oretically motivated) to why this is also a temporal 
regime inaccessible for exploration by traditional aver- 
aging schemes; expanding the temporal capabilities of 
single molecule studies. 

To bolster this position, we also presented mock (nu- 
merical simulations) experiments, demonstrating the fail- 
ure of traditional averaging methods to predict both 
the mechanical and kinetic properties of biomolecular 
systems, when the transitions between the biomolecu- 
lar states occur on the time scales governing the mo- 
tion of the experimental probe. To further motivate the 
attributes of the auto-correlation function, we also pre- 
sented numerical solutions which demonstrated how its 
functional form alters in this temporal regime. To com- 
prehend how the auto-correlation function may be influ- 
enced by the rates of transition, we then presented two 
physically motivated guesses for its solution. We argued 
that these two guesses should correspond to solutions 



for the auto-corrclation function in two physically oppos- 
ing extremes: when the switching between biomolecular 
states is either much faster or much slower than the time 
scales governing the motion of the probe. 

Motivated by these numerical simulations we then pre- 
sented our mathematical analysis for the solution to the 
auto-correlation function. We modeled the experimental 
probe as an over-damped particle, in which the associated 
spring constants corresponding to the biomoleculc fluc- 
tuates between different states. We presented a formal 
solution for the auto-correlation function, from which we 
proved that the two physically motivated solutions also 
correspond to mathematical bounds (opposing) to the 
general solution. 

Although no known "closed" form solution exists for 
the problem at hand, we presented series solutions for the 
auto-correlation function. Inspired by the bounds found 
for the general solution, two series solutions for the auto- 
correlation function were expressed around these bounds; 
providing series solutions which converge to the exact 
solutions from opposing positions. In addition, we proved 
mathematically that the two bounds indeed correspond 
to exact solutions for the auto-corrclation function in the 
anticipated regimes. 

Having obtained two series solutions for the auto- 
correlation function we next presented a general anal- 
ysis regarding their convergence to the general solution 
and their applicability to the analysis of single-molecule 
switching experiments. To garner concrete understand- 
ing regarding the convergence and applicability of this 
method, we then focuses on the two state problem. A 
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range of biomolccular parameters, for where we believe 
this method would be applicable for the analysis of ex- 
perimental data, was then determined. Guided by this 
analysis, we then presented solutions for each series (up 
to third order) for a two state Markov process. We com- 
pared our truncated series solutions with numerical solu- 
tions for the auto-correlation function and showed that 
these truncated series solutions converged to the numer- 
ical solution to the degree anticipated from our general 
analysis regarding convergence. These truncated series 
solutions for the auto-correlation function are amiable to 
the analysis of experimental data, in the regime where 
traditional averaging procedures fail. 

The general analysis presented here should be appli- 
cable to single-molecule experimental studies. Although 
the range of parameters found for applicability of this 
method appears to be limited, this parameter range is 
beyond the reach for analysis with current traditional 
methods, providing an expanded region of phase space 
that can be explored by single- molecule methods. 
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APPENDIX A: VARIANCE OF A TIME 
AVERAGED PROPERTY 



We now present an argument determining the reduc- 
tion of the width of the distribution function for time 
averaged values (a;)T. The widths of the distribution 
functions are related to the square root of the variance, 
cr^. Consider measurements of a property x over a given 
time window T and the average value of this measure- 
ment, (a;)T, is tabulated. The relationship between the 
data set corresponding to the time averaged value to that 
of the data set corresponding to the instantaneous mea- 
surements is characterized by the following equality 



data set), 
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(Al) 



Here Xi^j corresponds to an instantaneous measurement 
of X and the time averaged value corresponds to the 
average of x over time window i. Each time window 
corresponds to M data points and there are a total of 
NxM instantaneous measurements. 

Now we make a comparison between the variance cor- 
responding to the instantaneous measurements (whole 
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with that of the time averaged value, 
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Here x is the average value of x, 

N,M 



(A4) 



A few assumptions are now made. First it is assumed 
that the instantaneous measurements are uncorrelated, 
statistically independent. This implies that the instanta- 
neous measurements are stored at time intervals on order 
of the intrinsic time constant of the experimental probe, 
as this is roughly the time scale to make two uncorre- 
lated measurements. Next, we assume that the size of the 
data set corresponding to the time averaged value (N) is 
"large" such this data set {{{x').}) adequately samples 
the distribution space of the time averaged value. Fur- 
thermore, we assume that the number of data points used 
to compute each averaged value (M) is "small" such that 
the distribution function corresponding to the time av- 
eraged value contains a measurable width, cr^^,^^ ^ 0. 
(This also implies N^M.) Finally, we assume that all 
distribution functions have mirror symmetry around the 
mean (x) of the measured value x. These conditions are 
consistent with many single molecule experiments. 

Given these assumptions we are now in position to 
relate the variance corresponding to the instantaneous 
measurements, ct^, to that corresponding to the time av- 
eraged measurements, c^^^^ . The variance corresponding 



to the time averaged value, as expressed in Ea. lA3l can 
be related to the instantaneous measurements as follows. 




(A5) 
(A6) 

• (A7) 

Because of the assumption that M is small, implying that 



M 



j/M X, the last term does not vanish. This is 



consistent with having a nonzero variance corresponding 
to the time averaged value, ^(j;)^ ^ 0. 

Further decomposing, the variance fEa. lA7p can be 
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written as 



^ N M 

i j 
N M,M 



(A8) 



NM 



£ i^^'i - ^) (^''fc " ^) ■ (^9) 



Next, we use the assumption that N is "large". In the 
hmit that N goes to infinity the last term in Eq.EHl van- 
ishes, 



lim y](xij - x) ixi^k - a;) = 0, 



(AlO) 



for j ^ k. This equality holds when the data set {x)^ 
contain a number of statistically independent points and 
the distribution function has mirror symmetry around 
the average value x. 

The variance corresponding to the time averaged value 
now can be written as 
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Revealing that the variance corresponding to the time 
average value is equal to the variance corresponding to 
the instantaneous measurements, reduced by the num- 
ber of uncorrelated measurements used to determine the 
time averaged value. 



APPENDIX B: TIME CONSTANT 

The time scales governing the motion of the probe cor- 
respond to its time constant r = 7/fc, where k is the 
spring constant associated with the end of the probe. 
This time scale is also the time scale at which two un- 
correlated measurements are made. This spring constant 
contains contributions from both the intrinsic time con- 
stant of the experimental probe, fcp, and that from the 
molecule in state i, /cmi- This time constant is less than 
the intrinsic time constant of the experimental probe, 
'Tp ~ j/kp. For systems of interest the spring constant 
of the probe is generally much larger than or of order 
of the molecule in any state and; therefore, the intrin- 
sic time constant of the experimental probe serves as an 
appropriate measure for the relaxation time of the probe 
when attached to a biomolecule, rp « r. Using the in- 
trinsic time constant of the probe as the measure of its 
relaxation also proves to be a convenient measure, as the 



particular state of the molecule docs not have to be con- 
tinually referred to and quantified during this analysis. 
(Certainly, the discussions here can be generalized to the 
cases when the spring constant of the molecule in any 
state is much larger than the intrinsic spring constant of 
the probe.) 



APPENDIX C: AUTO-CORRELATION 
FUNCTION, WHEN (/(ti)/(t2) ■ ■ ■ /(tn)> IS 
INDEPENDENT OF TIME 



If all of the moments of the random variable / are 
independent of the time variables then Eq.221 takes the 
form 



{5x{t)Sx{0)) 
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where 



Gn+n' = (^„(0,{Tj}).F„.(t, {Tj,})) 



is independent of the time variables. It, therefore, only 
depends on the sum of n and ??' and not their individual 
values. 

To evaluate this integral wc first calculate the Fourier 
transform of the auto-correlation function and then in- 
tegrate over the time variables. We will then take the 
inverse transform by using the residue theorem for com- 
plex integration. Finally, through algebraic identities we 
obtain the result in Eg .1511 

The Fourier transform f{'jj) of a function f{t) is defined 



as 



dtf{t)e- 



and the inverse transform is defined as 

The Fourier transform of the auto-correlation function in 
this limiting form is 



1 \ '^Cn^n' I '^a 



n,7i'—0 \ " / 

Integrating over all the time variables, {ti} and {t;'}, the 
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Fourier transform becomes 



Using the identity 
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Taking the inverse transform of Ea. lC2l the auto- 
correlation function is 
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where Xl-'l is the j*'' derivative of X , and switching vari- 
ables to m = n + ?i' and rt' , Ea. lC4l can be written as 
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Using the identity 
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Equation IC5I can be written as 



Equation IC3I can be integrate in the upper-half 
complex-plane. The integrand has a pole at = ia/Ta 
and 
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the desired result. 



APPENDIX D: EXPANSION FOR A TWO STATE 
MARKOV PROCESS 

We present solutions for the series expansions repre- 
sented in Eas. l54l and 1551 up to third order, for a system 
that behaves as a two state Markov process. Here, a 



biomolecule is considered a two state Markov process if 
two distinct mechanical states exist (different spring con- 
stants) and the probability of making a transition from 
one state to another is local in time (independent of all 
previous configurations). We provide closed form solu- 
tions for the auto-correlation function (truncated at third 
order) that depends on the magnitude of the rates and 
can be used in the analysis of experimental data. We 
proceed by first reviewing the relevant theory of a two 
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state Markov process. 



1. The two state problem: Statistics 

We review some statistics governing a two state 
Markov process. Particularly, properties required for the 
evaluation of the moments {T>^, , 1^^°^/ ) expressed in the 
solution of the auto-correlation function, Eas.l54l and 1551 

Consider a two state Markov process and assume that 
at time i = the system starts in state i 0|. The 
probability of the system lying in state j at time t obeys 
the following rate equation 



dt \ P^2{t) 



-Ki2 H2l 
Kl2 — K21 



P,2(t) 



(Dl) 



Here Pij (i) is the probability of the system lying in state 
j at time t if at time t = Q the system was in state z, K21 
is the rate of going into state 1 from state 2 and K12 is 
the rate of going into state 2 from state 1. 
Define the rate matrix as 



W = 



^12 



1^21 



The probability of the system lying in a particular state 
{Pij ) at time t follows from solving the differential equa- 
tion represented in Ea. lDll In matrix notation the solu- 
tion is 



Pit) = TtP(O), 



(D2) 



Here the transition probability matrix, Tj , is defined as 



(D3) 



whose component, Tt^ij, is defined as the probability 
of the system lying in state i at time t if it was in state j 
at time 0. Because we assume a Markov process, the sta- 
tionary solution to Ea. lD2l and the transition probability 
matrix is all that is needed to obtain the statistics gov- 
erning the transitions of the system. All higher order, 
multi- variable, joint probability distributions functions 
can be obtained from them. 

The stationary solution to Ea. lD2l is 
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Pi 
P2 



(D4) 



where 



and 



Pi = K21/K, 
P2 = K12/K 



K = K.21+ K.12- 



It is straightforward to show that the transition proba- 
bility matrix is 



Tt = 1 - (1 - e-"*) 



P2 -Pi 
-P2 Pi 



where 1 is the identity matrix. 

For the moment we refine our discussion to instances of 
time. The stochastic variable of interest is the fluctuating 
term f{t), Eq.lJHl At time t„ the value of the fluctuating 
term is denoted as Z^"-*, which can take one of two values: 
/i if in state 1 and /2 if in state 2. For the problem at 
hand these values are fi = =Fl, where the minus sign cor- 
responds to the system with the larger spring constant, 
(see Ea. l36|) . Because we now consider instances of time 
we define slight variants to the stationary probability so- 
lution, P, and transition probability matrix, Tj. These 
variants refer to particular values of / instead of particu- 
lar states of the system. The stationary probability of / 
having value Z^"-* is denoted as Pi(/^"-') (single variable 
probability distribution function) . The probability of / 
having value /^"-' at time if it have had value /^'^■' at 
time tp is denoted as T*^^ (/(") l/f/^)). 

The multi-variable joint probability distribution func- 
tions, defined as the probability of / having value 
at time t„ if / had value /("~^) at time i„_i and value 
j(n-2) time tn-2 and so on, is 

P„(/(l\ti;/(2),i2;... ;/(«), = 

xPi(/(i)). 

Here time is ordered, <i < ^2 < < Ul, and tm.m-i = 
tm — trn-i > is the time difference of neighboring time 
instances. 

Averages over the stochastic variable / are needed to 
evaluate the integrals in Eqs.|^ and [S^ They can be 
written as 

{f{tl)f{t2)-.f{tn)) = 

/(l'-/(">={/l,/2} 

or in an alternative matrix form 

/n— 1 

{f{tl)f{t2)-.f{tn)) = (/l /2) n(Tl+e^''*-' 

where ti < t2 < ■■■ ^ tn (time ordering). Here 
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Given the averages of the fluctuating term, Ea. lD6l terms 
in the series expansion can be evaluated for the general 



27 



der. Expanding around the fast-switching hmit the auto- 
correlation function takes the form 



Truncated series 



Using the definitions above, we calculated the gen- 
eral solution for the auto-correlation, up to third or- 
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where the superscript ((00)) denotes that the scries is ex- 
panded around the fast-switching limit and the subscript 
denotes the order at which the series is truncated. In 
Ea. lDQI we present new notation. The vector r„(i) is a n- 
component vector that is a function of time and defined 
as 
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The coefficient vectors, C„ and I?„, are functions of kt^, 
whose values are listed in Table Ull Dagger (f ) represents 



the transpose of a vector and (/) = /i-Pi + /2-P2 is the 
average value of the fluctuating variable. The subscripts 
on the vectors T„(t), C„ and D„ represents the order 
each term was calculated, second or third, and also cor- 
responds to the size of the vector. In limit that kTq — > 00 
the coefficient vectors in Eg . ID9I vanish, leaving only the 
fast-switching limit solution. 



As expressed in Sect. once the fast-switching limit 
is determined expansion around the slow-switching limit 
is easily determined by exploiting Eas. lDQI and 1581 



where the superscript ((0)) denotes that the series is ex- 
panded around the slow-switching limit. Here, in the 
limit that kt^ the terms in the square brackets sum 
to zero, leaving the slow-switching limit solution. 

Equations ID9I and IDllI are general solutions for the 
auto-corrclation function. These solutions depend on the 



(62 + e-'=*i32 - 1) ^ f2it) + ^(/) (4 -I- e-'=*^3 - ^^T^{t) 

. Z mi 



magnitude of the rates of transition for the molecule, K12 
and K21, unlike the solutions in the fast-switching and 
slow-switching limits, which depend only on the ratio of 
the rates. 
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